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22 Abstract 

23 Bats harbor severe acute respiratory syndrome-like coronaviruses (SL-CoVs) from 

24 which the causative agent of the 2002-3 SARS pandemic is thought to have originated. 

25 However, despite a large number of genetically diverse SL-CoV sequences detected 

26 in bats, only two strains (named WIV1 and WIV16) have been successfully cultured 

27 in vitro. These two strains differ from SARS-CoV only in containing an extra ORF 

28 (named ORFX) between ORF6 and ORF7 with no homology to any known protein 

29 sequences. In this study, we constructed a full-length cDNA clone of SL-CoV WIV1 

30 (rWIV 1), an ORFX deletion mutant (rWIV 1 -AX), and a GFP-expressing mutant 

31 (rWIVl-GFP-AX). Northern blot and fluorescent microscopy indicate that ORFX was 

32 expressed during WIV1 infection. Virus infection assay showed that rWIVl-AX 

33 replicated as efficiently as rWIVl in Vero E6, Calu-3, and HeLa-hACE2 cells. 

34 Further study showed that ORFX could inhibit interferon production and activates 

35 NF-kB. Our results demonstrate for the first time that the unique ORFX in the WIV 1 

36 strain is a functional gene involving modulation of host immune response, but not 

37 essential for in vitro viral replication. 
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38 Importance 

39 Bats harbor genetically diverse SARS-like coronaviruses (SL-CoVs) and some of 

40 them have the potential of interspecies transmission. A unique open reading frame 

41 (ORFX) was identified in the genome of two recently isolated bat SL-CoV strains 

42 (WIV1 & 16). It will therefore be critical to clarify whether and how this protein 

43 would contribute to virulence during viral infection. Here we revealed that the unique 

44 ORFX is a functional gene involving in the modulation of host immune response, but 

45 not essential for in vitro viral replication. Our results provide important information 

46 for further exploration of the ORFX function in the future. Moreover, the reverse 

47 genetics system we constructed will be helpful for pathogenesis study of this group of 

48 viruses and develop therapeutics for future control of emerging SARS-like infections. 

49 

50 Introduction 

51 Severe acute respiratory syndrome coronavirus (SARS-CoV) is a zoonotic pathogen 

52 that caused the SARS pandemic in 2002-3 which originated in China (1). Since then, 

53 genetically diverse SARS-like coronaviruses (SL-CoVs) have been reported in bats in 

54 China, Europe and Africa (2-11), indicating a wide geographic distribution of this 

55 group of viruses. However, most bat SL-CoVs have only been identified by sequences 

56 and are not fully characterized due to the lack of cultured viruses. Thus, their potential 

57 for transmission to, and their likely pathogenesis in domestic animals and humans 

58 remain untested. WIV1 &16 are two recently identified SL-CoV strains with high 
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59 genomic similarity to human SARS-CoV. These two strains have been successfully 

60 cultured in vitro and have been shown to use the same molecule (angiotensin 

61 converting enzyme, ACE2) for cellular entry as SARS-CoV (2, 10). Recently, another 

62 bat SL-CoV strain SHC014 has been demonstrated to use human ACE2 by the 

63 construction of infectious cDNA clone (12). Furthermore, animal infection 

64 experiments indicated that SL-CoV WIV1 and SHCOM could replicate efficiently and 

65 caused low pathogenesis in ACE2 transgenic mice (12, 13). The fact that the native 

66 bat SL-CoVs could use human ACE2 without any mutations indicates a high risk of 

67 interspecies transmission for these and similar coronaviruses that may exist in natural 

68 reservoirs. 

69 Coronaviruses have the largest genomes among RNA viruses. Their genome 

70 consists of a positive, single-stranded RNA around 30 kilo base nucleotides (nt) with 

71 two-thirds at the 5’-end encoding genome replication proteins (ORFlab) and 

72 one-third at 3’-end encoding structural proteins including a spike glycoprotein (S), a 

73 small envelope protein (E), a membrane protein (M) and a nucleocapsid protein (N). 

74 Coronaviruses encode a set of ORFs expressed from full-length and 

75 subgenomic-length mRNAs (sgRNAs), which have common 3’ end originating at 

76 distinct transcription regulatory sequences (TRS) and joined with a common leader 

77 sequence encoded at the 5’ end of genomic RNA (14). Currently, coronaviruses are 

78 divided into the genera alphacoronavirus, betacoronavirus, gammacoronavirus, and a 

79 proposed genus deltacoronavirus (15). SARS-CoV and SL-CoVs are grouped into the 
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80 same coronavirus species, SARS-related coronavirus (SARSr-CoV), within the genus 

81 Betacoronavirus. Besides the family-conserved genes, SARSr-CoV possesses several 

82 accessory genes including ORF3, ORF6, ORF7, ORF8 and ORF9 which are specific 

83 for this group of coronaviruses, but not essential for in vitro viral replication (16-18). 

84 Accessory genes in coronavirus genomes play important roles in regulating host 

85 immune response (19). The SARS-CoV ORF3a, ORF3b, and ORF6 have been 

86 reported to inhibit host interferon (IFN) response during virus infection and contribute 

87 to pathogenesis (20, 21). ORF3a and ORF7a activate NF-kB and upregulates IL8 and 

88 CCL5 production (22, 23). Bat SL-CoVs display great genetic diversity and share 

89 overall nt sequence identities of 88-95% with human SARS-CoV (2-11). Bat SL-CoV 

90 WIV1 & 16 are the closest relatives, so-far discovered, to human SARS-CoV. These 

91 two viruses are identical in genomic structures except that WIV1 &16 have an extra 

92 ORF (named ORFX) between ORF6 and ORF7 with no homology to any known 

93 protein sequences (2). 

94 In this study, we explored the function of ORFX in modulating host immune 

95 response through eukaryotic overexpression assays and recombinant viruses generated 

96 through reverse genetics techniques. 

97 

98 Materials and Methods 

99 Virus and Cells. The SL-CoV WIV1 strain (GenBank accession: KF3 67457) and 
100 other viruses were propagated as described previously (2). Sendai virus (SeV) strain 
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101 Cantell (kindly provided by Prof. Hanzhong Wang) was propagated in 10-day-old 

102 embryonated chicken eggs at 37 °C for 48 h (24). All experiments using live virus 

103 was conducted under biosafety level (BSL) 2 conditions. HeLa cells stably expressing 

104 human ACE2 (HeLa-hACE2) were described previously (25). 293T, Vero E6, HeLa 

105 and HeLa-hACE2 cells were grown and propagated in Dulbecco’s Modified Eagle’s 

106 Medium supplemented (GIBCO, Invitrogen) with 10% fetal bovine serum (Life 

107 Technologies). Calu-3 cells were grown and propagated in Dulbecco's Modified Eagle 

108 Medium: Nutrient Mixture F-12 Media supplemented with 15% fetal bovine serum. 

109 Cells were grown at 37 °C in a humidified atmosphere with 5% CCK 

110 

111 Plasmids. The encoding region of the ORFX was amplified by reverse transcription 

112 (RT)-PCR from viral RNA using Superscript One-Step RT-PCR kit (Invitrogen). The 

113 amplified gene was cloned into the plasmid pCAGGS with a C-terminal HA tag 

114 (pCAGGS-ORFX) for eukaryotic expression. Reporter plasmids used included 

115 pIFN(3-Luc (expressing firefly luciferase under the control of the IFN-(3 promoter), 

116 pNF-KB-Luc (expressing firefly luciferase under the control of the NF-kB promoter) 

117 and pRL-TK (expressing Renilla luciferase under the control of the HSV-TK 

118 promoter), as well as expression plasmid for influenza virus NS1 as described 

119 previously (24). Subcellular organelle markers expressing plasmids including 

120 SecGip-GFP (Endoplasmic Reticulum (ER) marker), B4Gal-Ti-RFP (Golgi marker), 
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121 and Mito-YFP (mitochondria marker) were kindly provided by Prof. Yanyi Wang of 

122 the Wuhan Institute of Virology. 

123 

124 Viral infection assays. Vero E6, Calu-3, and HeLa-hACE2 cells were infected with 

125 viruses at a multiplicity of infection (MOI) of 1.0, 0.1, or 0.001 in 25-cm 2 flasks with 

126 1 h adsorption period, followed by two washes with D-Hanks and culturing by adding 

127 3 mL of medium. The viral supernatants were harvested, at 0, 2, 6, 12, 18, 24, 36, 48, 

128 and 72 h post inoculation, with 300 pL removed and 300 pL medium added back at 

129 each time point. Virus concentration was titrated by plaque assay in Vero E6 cells. 

130 Vero E6 cells were infected by rWIVl-GFP-AX or mock infected. After 24 h, 

131 fluorescent micrographs was taken to check the expression of green fluorescence 

132 protein. 

133 

134 Cloning of WIV1 cDNAs. The virus genome was divided into 8 continuous 

135 fragments (A-G) and amplified using specific primers (primer sequences available 

136 upon request). Viral RNA was extracted from the supernatant of WIVI-infected 

137 cultures, and reverse-transcribed with M-MLV reverse transcriptase (Promega), 

138 random hexamer deoxynucleotide primers. The cDNA was denatured for 5 min at 

139 95 °C and amplified by polymerase chain reaction (PCR) with KOD DNA 

140 polymerase (TOYOBO) for 20 cycles of 95 °C for 30 sec, 60 °C for 30 sec with 

141 0.5 °C decrease per cycle, 68 °C for 5 min, and 15 cycles of 95 °C for 30 sec, 50 °C 
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142 for 30 sec, 68 °C for 5 min, and a final extension at 68 °C for 10 min. The amplicons 

143 were cloned into pGEM-T EASY (Promega). Besides three natural Bgll sites, several 

144 Bgll sites were introduced by synonymous mutations in the PCR process, to make all 

145 contiguous cDNA fragments capable of unidirectional ligation. SacII and AscI were 

146 introduced into the 5’ terminus of fragment A and the 3’ terminus of fragment G 

147 respectively. A poly(A) sequence (25 nt) was added to the 3’ terminus of fragment G. 

148 At least three colonies of each cDNA clone were sequenced, and the one identical to 

149 or with some synonymous mutations to the reported sequence was selected for 

150 assembly. 

151 To ablate a natural Bgll site at position 1575, primers FA, F-cl575a, R-cl575a 

152 and RA were used for overlap extension PCR (OE-PCR) to introduce a synonymous 

153 mutation C1575A (primer sequences available upon request). Based on previous in 

154 vitro transcription tests, a synonymous mutation T27527C was also introduced to 

155 interrupt a potential T7 termination site via OE-PCR. 

156 

157 Strategy for modifying pBeloBACl 1. 

158 Cytomegalovirus (CMV) promoter was amplified from pcDNA3.1(+) (Thermo 

159 Fisher Scientific) with primers (forward, 5’-tgaggatcccgttgacattgattattgactag-3’; 

160 reverse, 5’-cctgactgcaggtcgactgccgcggagctctgcttatatagacc-3’). Flepatitis delta virus 

161 (F1DV) ribozyme was synthesized as described (26), and amplified with primers 

162 (forward, 5’-cagtcgacctgcagtcaggcgcgccgggtcggcatggcatctcc-3’; reverse, 
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163 5’-ctagaaggcacagctcccttagccatccgagtgg-3’). Bovine growth hormone (BGH) 

164 transcription terminal signal was amplified from pcDNA3.1(+) with primers (forward, 

165 5’-ggatggctaagggagctgtgccttctagttgccagc-3’, reverse, 

166 5’-tgaaagcttccatagagcccaccgcatcc-3’). Then the three PCR products were ligated 

167 using OE-PCR, with BamHI and Hindlll sites flanking the amplicon, and SacII and 

168 AscI sites sitting between CMV promoter and HDV ribozyme. The amplicon was 

169 then inserted into pBeloBACl 1 (New England BioLabs) between BamElI and Hindlll 

170 sites. The construct was designated pBAC-CMV. 

171 

172 Construction of infectious BAC clones of WIV1. The subclone A and subclone G 

173 were first digested with SacII and AscI (New England BioLabs) respectively, 

174 followed by treatment with CIAP (Takara), chloroform extraction and isopropanol 

175 precipitation, and then restricted with Bgll (Takara). The subclones B-F were digested 

176 with Bgll. pBAC-CMV was digested with SacII and AscI. All digestion products 

177 were then separated using 1% agarose gels, excised, and purified by using the Gel 

178 Extraction Kit (Omega). Digested fragments A-G and pBAC-CMV were ligated 

179 overnight at 4 °C, transformed into DH10B competent cells and plated on Chl + LB 

180 culture. Ten clones were screened by RFLP analysis with Ncol, StuI or Hindlll. The 

181 correct clone was named as pBAC-CMV-rWIVl (Fig. 1). 

182 
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183 Construction of WIV1 mutants. To delete ORFX, the fragment F was 

184 PCR-amplified with primer set FF (5’-acctgtgcccttttggcgaggtttttaatgctactac-3’) and 

185 RFox (5’-gcctctagggctcaaggataatctatctccatagg-3’). The fragment G was 

186 PCR-amplified with FGox (5’-gccctagaggcaacgaacatgaaaattattctcttcc-3’) and RG (5’- 

187 actggcgcgcctttttttttttttttttttttttttgtcattctcctgagaagc-3 ’). This new fragment is named as 

188 Gox. These two products were then cloned into pGEM-T EASY. The two fragments 

189 were inserted into BAC along with the other fragments as described above. The 

190 rescued mutant was named as rWIVl-AX. To replace GFP into the open reading 

191 frame of ORFX, the F fragment was amplified with primer set FF and RFoeGFP 

192 (5’-gctcaccatagtggttcgtttatcaaggataatctatctcc-3’). The GFP gene was amplified with a 

193 primer set (5’-ccttgataaacgaaccactatggtgagcaagggcgaggagc-3’ and 

194 5’-tgcctctagggcttacttgtacagctcgtccatgcc-3’). The two PCR products were ligated by 

195 OE-PCR and the product was inserted into pGEM-T EASY. The rescued mutant was 

196 named as rWIVl-GFP-AX. 

197 

198 Transfection of infectious WIV1 BAC clones. Vero E6 cells were seeded in a 6-well 

199 plate a day in advance, and then one well was transfected with 6 pg infectious BAC 

200 plasmids constructed as above with Lipofectamine LTX and Plus Reagent (Life 

201 technologies). Virus progeny was plaque purified once. One clone was passaged once 

202 in Vero E6 cells for 72 h and used to generate a stock for future use. 

203 
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204 Restriction fragment length polymorphism. RNA extracted from wild-type and 

205 recombinant viruses were reverse-transcribed with random hexamer primers. RT-PCR 

206 was used to generate five amplicons containing the five mutations designed in the 

207 strategy. These amplicons included a 1124-bp amplicon (nucleotide positions 

208 1312-2435) spanning a naturally occurring Bgll site at nucleotide 1571 that had been 

209 ablated in recombinant viruses, a 1438-bp amplicon spanning the B/Cl junction 

210 (nucleotide positions 7560-8997), a 1437-bp amplicon spanning the C1/C2 junction 

211 (nucleotide positions 10196-11632), a 1437-bp amplicon spanning the D/E junction 

212 (nucleotide positions 16793-18229), and a 1438-bp amplicon spanning the E/F 

213 junction (nucleotide positions 21908-23345) (These amplicons correspond to 

214 fragment F1-F5 in Fig. 1). The first amplicon of wild-type WIV1 that contains 

215 nucleotide 1571 can be cleaved by Bgll, but the other four amplicons cannot. In 

216 contrast, the five amplicons of recombinant viruses are different to those of wild-type 

217 virus in the capability of being cut by Bgll. 

218 

219 Northern blot analysis. The N gene was amplified with primers WIV1-NF 

220 (5’-atgtctgataatggacccca-3’) and WIV1-3R (5’-gtcattctcctgagaagcta-3’) and used as a 

221 template for probe preparation according to the description of DIG-Fligh Prime DNA 

222 Labeling and Detection Starter Kit II (Roche). Vero E6 cells were infected with 

223 wild-type and recombinant viruses at an MOI of 1.0. At 24 h postinfection, 

224 intracellular RNA was isolated using TRIzol reagent (Ambion). RNA (20 pg) was 
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225 precipitated, treated with 17 pL Sample Buffer (50% Formamide, 2.2 M 

226 Formaldehyde (37%), lx MOPS) at 65 °C for 10 min, added with 3 pL 10x Dye 

227 Solution (50% glycerol, 0.25% Bromophenol Blue, 0.25% Xylene Cyanole FF), and 

228 then separated in denaturing 0.8% agarose-2.2 M formaldehyde gel at 28 V for -17 h. 

229 The RNA was hydrolyzed by 0.05 M NaOFl for 40 min, transferred to a Flybond-N+ 

230 membrane (GE Flealthcare) for ~18 h, and then cross-linked to the membrane using 

231 UV light. The membrane was prehybridized, probed with a DIG-labelled probe for N 

232 gene, washed and detected according to the DIG-Fligh Prime DNA Labeling and 

233 Detection Starter Kit II (Roche). 

234 

235 RT-PCR of leader-containing transcripts. Intracellular RNA was isolated from 

236 wtWIVl. A forward primer (Leader-F) located in the leader sequence, along with 

237 various reverse primers located in several ORFs, was used for amplifying 

238 leader-containing sequences (primer sequences available upon request). 

239 Leader-containing amplicons were sequenced with the corresponding reverse primers. 

240 

241 ORFX subcellular location. FleLa cells were transfected with ORFX expressing 

242 plasmid and cotransfected with organelle markers expressing plasmids SecGl|3-GFP, 

243 B4Gal-Ti-RFP, or Mito-YFP. After 24 h, the cells were fixed and stained with a 

244 mouse anti-HA IgG (Promoter). A Cy3 conjugated goat anti-mouse IgG (Promoter) 

245 was used for secondary detection in cells expressing ER or mitochondrial markers. An 
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246 FITC conjugated goat anti-mouse IgG (Promoter) was used for secondary detection in 

247 cells expressing Golgi marker. Nuclei were stained with 

248 4',6-diamidino-2-phenylindole (DAPI). Staining patterns were examined with an 

249 Olympus Fluoview upright confocal microscope (Olympus). 

250 

251 Luciferase assays and quantitative PCR. For ORFX-mediated IFN promoter assay, 

252 293T cells were seeded in 12-well plates and cotransfected with empty vector plasmid 

253 pCAGGS, plasmid pCAGGS-NSl or escalating amount (100, 200, 400, 600, 800 ng) 

254 of pCAGGS-ORFX with the indicated reporter plasmids. At 24 h posttransfection, 

255 cells were infected with SeV (100 hemagglutinin units/ml) for 12 h to induce IFN 

256 production, or treated with TNFa for 1 h to activate NF-kB. Cell lysates were 

257 prepared and luciferase activity was measured using the Dual-Luciferase Assay kits 

258 (Promega) according to the manufacturer’s instructions. 

259 293T cells were transfected with empty vector, NS 1 expressing plasmid, or 

260 escalating amount (100, 300, 600 ng) of ORFX expressing plasmid. After 24 h, the 

261 cells were infected with SeV (100 HAU/ml). At 12 h postinfection, the cells were 

262 lysed. The mRNA was extracted and reverse transcribed with PrimeScript RT Master 

263 Mix (Takara). The expression level of IFN-J3 mRNA was determined by 

264 quantification PCR using SYBR Premix Ex Taq II (Takara). The GAPDH mRNA was 

265 quantified as an inner control. 293T cells were transfected as above. After 24 h, the 

266 cells were treated with TNFa for 6 h, the cell RNA was extracted and used for 
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267 quantification of the expression of IL8 mRNA. All experiments were performed in 

268 triplicate and repeated at least three times. All primer sequences used in the 

269 quantitative PCRs will be provided upon request. 

270 

271 IRF3 translocation assay. 293T cells were transfected with empty vector, NS1, or 

272 ORFX expressing plasmid. After 24 h, IRF3 nuclear translocation was induced by 

273 infecting the cells with SeV for 8 h. The cells were fixed and stained with a rabbit 

274 anti-IRF3 polyclonal IgG (Proteintech) and a mouse anti-HA IgG (Promoter). An 

275 Alexa Fluor 488 conjugated donkey anti-rabbit IgG (Yeasen) and an Alexa Fluor 555 

276 conjugated donkey anti-mouse IgG (Beyotime) were used to detect IRF3 and ORFX, 

277 respectively. The cells transfected with empty vector were stained with a rabbit 

278 anti-IRF3 polyclonal IgG and a goat anti-SeV IgG (kindly provided by Prof. Lin-Fa 

279 Wang at the Duke-NUS Graduate Medical School, Singapore), as an indication of 

280 infection efficiency. An Alexa Fluor 488 conjugated donkey anti-rabbit IgG and a 

281 Cy3 conjugated donkey anti-Goat IgG (Promoter) were used to detect IRF3 and SeV, 

282 respectively. Nuclei were stained with DAPI. 

283 

284 Quantification of mRNA expression of cytokines in infected Calu-3 cells. Calu-3 

285 cells grown in 24-well plates were mock infected, or infected with rWIVl or 

286 rWIVl-AX at an MOI of 5, or SeV (100 FLAU/ml). The cells were lysed at 4, 12, 24, 

287 and 30 h postinfection. The mRNA expression level of IFN-J3, IL6, IL8, and TNFa 
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288 was quantified by quantitative PCRs. The expression of GAPDH mRNA was 

289 measured as an inner control. All primer sequences used in the quantitative PCRs will 

290 be provided upon request. The experiment was performed twice. 

291 

292 IFN-P sensitivity assay. Vero E6 cells were seeded a day in advance. The cells were 

293 pretreated with 10, 100, or 1000 U/ml IFN-P (PBL, Piscataway, NJ) for 24 h, infected 

294 with wtWIV 1, rWIV 1, and rWIV 1 -AX at an MOI of 0.1 PFU/cell, and post treated 

295 with the same amount of IFN-P as used previously. At 24 h postinfection, the viral 

296 replication was analyzed by plaque assay. The experiment was performed in triplicate. 

297 

298 Statistics. The statistical significance of the obtained data was analyzed using a 

299 student’s t-test in GraphPad Prism (GraphPad Software, San Diego, CA). A p-value 

300 <0.05 was considered statistically significant. Data is presented as the means ± SEM. 

301 

302 Results 

303 Strategy for construction of an infectious WIV1 BAC. Originally, the genome was 

304 split into seven contiguous cDNAs (A~G) (Fig. 1 A and C). Due to the plasmid 

305 instability, fragment C was separated into two segments (Cl and C2). Besides three 

306 naturally-occurring Bgll (GCCNNNNJ.NGGC) sites, four Bgll sites were successfully 

307 introduced by synonymous mutations in the genome (Fig. 1 B). Different asymmetric 

308 3-nt overhangs at the junctions of each two contiguous fragments were created by 
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309 these Bgll sites. The eight fragments were then linked in one direction. A SacII site 

310 was added to the 5’ terminus of fragment A. A poly(A) sequence (25 nt) and an AscI 

311 site were added to the 3’ terminus of fragment G. A naturally occurring Bgll site at 

312 nucleotide 1571 was removed by synonymous mutations C1575A (Fig. IB). Other 

313 unexpected synonymous mutations also occurred including T1422C, T12984C, 

314 T14213C, T17130C, C17934T and T26068G. 

315 The plasmid pBAC-CMV was constructed by inserting with the cytomegalovirus 

316 (CMV) promoter, hepatitis delta virus (HDV) ribozyme and bovine growth hormone 

317 (BGH) transcription terminal signal sequences into pBeloBACl 1, along with the 

318 introduction of the SacII and AscI sites between CMV promoter and HDV ribozyme 

319 (Fig. 1 C). The eight genomic fragments were inserted into pBAC-CMV in one step. 

320 Recombinant viruses could be rescued by direct transfection with the BAC constructs. 

321 

322 Rescue of recombinant viruses. To rescue recombinant WIV1 (rWIVl ), fragments 

323 A and G were digested with SacII and AscI, respectively. Following calf intestinal 

324 alkaline phosphatase (CIAP) dephosphorylation, the two fragments, along with 

325 fragments B~F were digested using Bgll and inserted into pBAC-CMV between SacII 

326 and AscI sites in one step. The constructed clone (pBAC-CMV-rWIVl) was 

327 transfected into Vero E6 cells. Cytopathic effect was observed at 72 h posttransfection. 

328 The one ablated natural Bgll site and four introduced Bgll sites in the rescued viral 

329 genome were confirmed by restriction fragment length polymorphism (RFLP) 
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330 analysis with Bgll digestion (Fig. 2 A). Using this method, we also rescued an ORFX 

331 deletion mutant virus (rWIVl-AX) (Fig. 2 B lane 2), and a mutant with GFP replaced 

332 into the coding region of ORFX (rWIVl-GFP-AX) (Fig. 3 A). 

333 

334 ORFX is a functional gene not essential for virus replication. The one-step growth 

335 curve of the two rescued recombinant viruses (rWIVl-AX and rWIVl) and wild-type 

336 WIV1 (wtWIVl) determined by plaque assay showed that rWIVl-AX and rWIVl 

337 both replicated to the titers close to wild-type virus (Fig. 2 C). The expected set of 

338 appropriately sized 10 sgRNAs including sgRNA7 (ORFX) were observed in 

339 Northern blot analysis in cells infected with wtWIVl and rWIVl (Fig.2B lane 1 and 

340 lane 3). As expected, sgRNA7 was not observed in rWIVl-AX infected cells (Fig. 2 

341 B lane 2). Analysis of leader containing sequences indicated that all 10 sgRNA in 

342 wtWIVl share an identical core sequence ACGAAC (Table 1), which further 

343 confirmed that ORFX is expressed as sgRNA7. The GFP was expressed in 

344 rWIVl-GFP-AX infected cells further confirmed that the open reading frame of 

345 ORFX could be expressed (Fig. 3 A). Subcellular location analyses showed that the 

346 ORFX protein colocalized with the ER marker, but not with the Golgi and 

347 Mitochondria markers (Fig. 3 B). 

348 

349 ORFX protein inhibits the production of IFN-p. To determine whether ORFX 

350 inhibits the induction of IFN, 293T cells were transfected with plasmids pIFNp-Luc 
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351 and pRL-TK and the plasmid expressing ORFX, influenza virus strain PR8 NS1 

352 (positive control) or empty vector (negative control). As expected, SeV activated IFN 

353 production in cells transfected with empty vector. The positive control, influenza 

354 virus NS1 protein, dramatically inhibited the expression from the IFN promoter. 

355 ORFX protein exhibited inhibition effect, but the effect decreased while more ORFX 

356 protein was expressed (Fig. 4 A). Similar results were observed for IFN-J3 mRNA 

357 quantification (Fig. 4 B and C). 

358 IRF3 nuclear translocation assay was performed to see whether ORFX protein 

359 inhibits IFN production through inhibiting this process. 293T cells were transfected 

360 with empty vector, NS1, or ORFX expressing plasmid. After 24 h, IRF3 nuclear 

361 translocation was induced by infection of SeV for 8 h. The relative IRF3 translocation 

362 ratios were calculated for each group by counting the number of the IRF3 nuclear 

363 translocation cells (randomly selected at least 4 fields) divided by the number of total 

364 infected or transfected cells. The IRF3 nuclear translocation efficiency of each group 

365 was expressed as the percentage of their relative IRF3 translocation ratios to that of 

366 control (cells transfected with empty vector). As expected, NS1 strongly inhibited 

367 translocation of IRF3. While ORFX protein also showed inhibition of IRF3 

368 translocation, but less efficiently (Fig. 4 D and E). 

369 To further investigate the IFN inhibition activity of ORFX, the deletion mutant 

370 and wild-type recombinant virus were used to infect Calu-3 cells at an MOI of 5. 

371 Mock-infected cells were used as negative control. Calu-3 cells infected with SeV 
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372 were used as positive control. Samples were collected at 4, 12, 24, and 30 h 

373 postinfection. The relative expression of IFN-P mRNA was determined by 

374 quantification PCR and normalized to the expression of GAPDH mRNA. Compared 

375 to SeV, WIV1 recombinants induced low levels of IFN-P mRNA in Calu-3 cells (Fig. 

376 4 F). The ORFX deletion mutant induced a significantly higher level of IFN-|3 mRNA 

377 than wild-type recombinant virus in infected cells at 12 h postinfection, but no 

378 significant differences at 24, 30 h postinfection (Fig. 4 F). These results indicate that 

379 ORFX protein may play a role in antagonizing IFN only in early time during WIV1 

380 infection. 

381 

382 ORFX deletion mutant shows increased sensitivity to IFN-p. To further investigate 

383 the effect of ORFX on the viral sensitivity of IFN, we tested the replication 

384 efficiencies of wtWIV 1, rWIV 1 and rWIV 1 -AX in Vero E6 cells which were 

385 pretreated and posttreated with IFN-p. The replication of rWIVl-AX was evidently 

386 inhibited and reduced ~0.5 logs compared to wtWIV 1 and rWIVl at a concentration 

387 of 10 and 100 U/ml IFN-P (Fig. 4 G). Whereas, at a higher IFN-P concentration (1000 

388 U/ml), the titers of rWIVl-AX didn’t show an obvious decrease compared to wild 

389 type virus. We expected the ORFX deleted mutant would replicate less efficiently 

390 than the wild-type virus in IFN competent cells. Flowever, we did not find a 

391 significant difference when we grew the two viruses in Calu-3 and HeLa-hACE2 cells, 

392 even at a very low MOI of 0.001 (Fig. 5). 
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393 

394 ORFX protein activates NF-kB. NF-kB plays an important role in regulating 

395 immune response to viral infection and is also a key factor frequently targeted by 

396 viruses for taking over the host cell (27). Several proteins (Nspl, N, and 3a) encoded 

397 by SARS-CoV have both activities in IFN antagonizing and NF-kB activation (28). In 

398 this study, we also tested whether ORFX protein could activate NF-kB. 293T cells 

399 were transfected with pNF-KB-Luc, pRL-TK, and empty vector, NS1, or escalating 

400 amount (200, 400, 600 ng) of ORFX expressing plasmid. After 24 h, the cells were 

401 mock treated, or treated with TNFa for 6 h, and luciferase activity was determined. 

402 ORFX protein obviously activates NF-kB, no matter whether the cells were treated 

403 with TNFa or not (Fig. 6 A). Whereas, IL8 was upregulated only when the cells were 

404 treated with TNFa (Fig. 6 B). However, no significant difference was observed for 

405 IL6 and IL8 transcription levels between the rWIVl-AX and rWIVl infected Calu-3 

406 cells (Fig. 6 C and D). A significant difference was only observed for the induction 

407 of TNFa mRNA at the late time of virus infection in which the ORFX deletion mutant 

408 induced apparently less amount of TNFa mRNA (Fig. 6 E). 

409 

410 Discussion 

411 In this study, we have developed a fast and cost-effective method for reverse 

412 genetics of coronaviruses by combining two approaches developed by others (29, 30). 

413 Our method allows the genomes of coronaviruses to be split into multiple fragments 
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414 and inserted into a BAC plasmid with a single step. Recombinant viruses can then be 

415 efficiently rescued by direct transfection of the BAC constructs. As the genomes can 

416 be divided into multiple short fragments, mutations can be introduced into individual 

417 fragments easily (31). Using this method, we successfully rescued three recombinant 

418 viruses derived from SL-CoV WIV1 (rWIVl, rWIVl-AX and rWIVl-GFP-AX). The 

419 recombinant rWIVl and rWIVl-AX replicated to titers close to wtWIVl in Vero E6 

420 cells (Fig. 2 C), suggesting that the deletion of ORFX did not affect WIV1 replication 

421 in vitro. Northern-blot and fluorescent microscopy further confirmed that ORFX is 

422 transcribed as sgRNA7 and translated in virus-infected cells. These results 

423 demonstrated that unique ORFX in SL-CoV WIV1 is a functional gene but not 

424 essential for virus replication. We propose that the ORFX sgRNA is the template for 

425 the translation of a novel 11-kDa accessory protein of WIV1, bringing the total 

426 number of group-specific accessory proteins to ten. 

427 In previous studies, it has been proved that SARS-CoV group-specific accessory 

428 genes ORF3b and ORF6 inhibit host IFN production and/or signaling during virus 

429 infection and contribute to viral pathogenesis (20). It’s interesting to know if the 

430 ORFX has a similar function in IFN antagonizing. In this study, ORFX protein 

431 showed inhibition effect on IFN production, but the effect decreased when more 

432 ORFX protein was expressed (Fig. 4 A and B). Moreover, ORFX deletion mutant had 

433 significant lower inhibition effect on IFN production than wild-type recombinant 

434 virus in infected Calu-3 cells, but only at an early time after infection (Fig. 4 F). 
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435 Furthermore, the IFN sensitivity assay indicated that the ORFX deletion mutant was 

436 more sensitive to IFN-(3 (Fig. 4 G), suggesting that ORFX protein may participate in 

437 subverting the antiviral state stimulated by IFN-J3. All these results suggested that 

438 ORFX participates in the modulation of IFN response. Previous studies showed that 

439 the SARS-CoV ORF3a and ORF7a activate NF-kB and upregulate IL8 and CCL5 

440 production (22, 23). In our study, we also found that overexpressed ORFX can 

441 activate NF-kB through dual luciferase assay (Fig. 6 A). Furthermore, the level of 

442 TNFa mRNA induced by wild-type recombinant virus was significantly higher than 

443 those induced by ORFX deletion mutant, but only at the late stage of infection (Fig. 6 

444 E). These results indicated that ORFX also participates in activation of NF-kB. We 

445 noted that the IFN inhibition activity of ORFX was not dose-dependent and decreased 

446 when there was more ORFX expression. One possible hypothesis is that ORFX 

447 inhibits the IFN only at the early stage of infection. At the late stage, it activates 

448 NF-kB which in turn stimulate IFN expression, this leads to the attenuation of its IFN 

449 antagonist activity. 

450 Coronavirus was previously shown to induce unfolded-protein response (UPR) 

451 and ER stress in infected culture cells (32). Normally, ER is an active organelle for 

452 protein folding and modification. Loss of protein folding homeostasis would cause ER 

453 stress and induce UPR, leading to the activation of three ER stress transducers. These 

454 transducers work in a concert to attenuate translation and improve ER folding 

455 capacity to restore ER homeostasis (33). In this process, NF-kB is activated, and 
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456 apoptosis will be induced if ER stress prolong (32, 33). In this study, we observed that 

457 the overexpression of ORFX protein lead to cell death and the decrease of Renilla 

458 values (data not shown). This may imply that ORFX has cytotoxic effect and an 

459 influence on overall protein translation. We also found that ORFX colocalizes with 

460 ER marker. We hypothesize that ORFX may induce UPR and cause ER stress which 

461 would activate NF-kB and induce apoptosis, promoting viral release at late stage of 

462 infection. 

463 It should be noted that the IFN and NF-kB detection systems used in this study 

464 were derived from and performed in human cells. Since the innate immune system of 

465 bats is special and probably deficient in some aspects compared to human (34). Thus, 

466 it will be interesting to conduct the same studies in bat cells to determine whether 

467 ORFX protein will have the same profiles as those observed in the human cell system. 

468 The development of different Rhinolophus bat cell lines, which is the reservoir host of 

469 SL-CoV, will facilitate this research in the future. 

470 
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598 Figure Legends 

599 Fig. 1 Strategy for construction of an infectious WIV1 BAC clone. (A) Genomic 

600 structure of WIV1. (B) The mutations are indicated under the stars (★). Cl575A was 

601 used to ablate a natural Bgll site at nucleotide 1571 (V) t and T27527C was used to 

602 disrupt a potential T7 stop site. The others were for introducing Bgll sites ( T). (C) 

603 The WIV1 genome was split into eight contiguous cDNAs (A-G): A (nt 1-4387), B 

604 (nt 4388-8032), Cl (nt 8033-10561), C2 (nt 10562-12079), D (nt 12080-17017), E (nt 

605 17018-22468), F (nt 22469-27352), G (nt 27353-30309). Unique Bgll sites were 

606 introduced into the fragments by synonymous mutations to make these fragments 

607 capable of unidirectional ligation along with native Bgll sites in the genome. The 

608 original nucleotides were shown upside of the flanking sequences of corresponding 

609 fragments. A poly(A) sequence was added to the 3’ terminus of fragment G. CMV 

610 promoter, HDV ribozyme and BGH transcriptional terminal signal were inserted into 

611 pBeloBAC 11 between BamHI and Hindlll. SacII and AscI were introduced between 

612 CMV promoter and ribozyme. The fragments A-G were inserted into the pBAC-CMV 

613 plasmid in a single step. 

614 

615 Fig. 2 Recovery and characterization of recombinant viruses. (A) Restriction 

616 fragment length polymorphism. Amplicons flanking five mutated sites of wild-type 

617 and recombinant viruses were digested by Bgll. The first amplicon (FI) of wild-type 

618 virus can be digested by Bgll and its other four amplicons (F2~5) cannot. On the 
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619 contrary for amplicons of rWIVl, the first amplicon (FI) cannot be digested by Bgll 

620 and its other four amplicons (F2~5) can. M: DL2000 DNA ladder (Takara). (B) 

621 Detection of viral genomic transcription and replication by Northern blotting. Vero E6 

622 cells were infected with wild-type or recombinant viruses and intracellular RNA was 

623 extracted for Northern blot analysis. Lane 1, wtWIV 1; lane 2, rWIV 1 -AX; lane 3, 

624 rWIVl; lane 4, uninfected control. (C) Growth kinetics of wild-type and recombinant 

625 viruses. Vero E6 cells were infected with wtWIVl (■), rWIVl (O), or rWIVl-AX 

626 (A) at an MOI of 1.0 or 0.1 PFU/cell. Cell supernatants were taken at indicated time 

627 points postinfection and virus titers were determined by plaque assay in Vero E6 cells. 

628 

629 Fig. 3 Expression and subcellular location of ORFX protein. (A) The open reading 

630 frame of ORFX was replaced by GFP and the recombinant virus was rescued. Vero 

631 E6 cells were infected with the recombinant virus or mock infected. Green 

632 fluorescence was visualized at 24 h postinfection. (B) ORFX protein with an HA tag 

633 at C terminal was expressed in HeLa cells, along with SecGl|3-GFP (ER marker), 

634 Mito-YFP (mitochondria marker), or B4Gal-Ti-RFP (Golgi marker), respectively. 

635 The cells were fixed after 24 h and stained with a mouse anti-HA IgG. A Cy3 

636 conjugated goat anti-mouse IgG was used for secondary detection in cells expressing 

637 ER or mitochondria marker. An FITC conjugated goat anti-mouse IgG was used for 

638 secondary detection in cells expressing Golgi marker. ORFX protein showed a 

639 cytoplasmic distribution and colocalized with ER maker SecGlfl. 
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640 

641 Fig. 4 ORFX protein inhibits the production of type I interferon. (A and B) 293T 

642 cells seeded in 12-well plates were transfected with 100 ng pIFN-|3-Luc, 5 ng pRL-TK, 

643 and empty vector, influenza A NS1 expressing plasmid, or escalating doses (100, 200, 

644 400, 600, 800 ng) of ORFX expressing plasmid. Empty vector was added 

645 appropriately to ensure that cells in each well were transfected with the same amount 

646 of plasmids. The cells were infected with Sendai virus (100 hemagglutinating units/ml) 

647 at 24 h posttransfection. Samples were collected at 12 h postinfection followed by 

648 dual-luciferase assay. The results were expressed as firefly luciferase value 

649 normalized to that of Renilla luciferase. The relative expression of IFN-(3 mRNA was 

650 determined by quantitative RT-PCR and normalized to the expression level of 

651 GAPDF1 mRNA. (C) The expression of NS1 and ORFX proteins were analysis by 

652 western blot with an antibody against FlA-tag. The experiments were replicated for 

653 three times. (D and E) For IRF3 translocation assay, 293T cells were transfected with 

654 empty vector, NS1, or ORFX expressing plasmid. After 24 h, the cells were infected 

655 with Sendai virus to induce IRF3 nuclear translocation. The cells were fixed at 8 h 

656 postinfection and stained with anti-HA IgG. A Goat anti-Sendai virus polyclonal IgG 

657 was used to stain the cells transfected with empty vector. A rabbit anti-IRF3 

658 polyclonal IgG was used to label IRF3. The white arrow indicates IRF3 nuclear 

659 translocation. The relative IRF3 translocation ratios were calculated for each group by 

660 counting the number of the IRF3 nuclear translocation cells (randomly selected at 
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661 least 4 fields) divided by the number of total infected or transfected cells. The IRF3 

662 nuclear translocation efficiency of each group was expressed as the percentage of 

663 their relative IRF3 translocation ratios to that of control (cells transfected with empty 

664 vector). (F) Calu-3 cells were mock infected, or infected with rWIVl, rWIVl-AX 

665 (MOI of 5), or SeV (100 HAU/ml). At 4, 12, 24, and 30 h postinfection, the cell RNA 

666 was extracted and used for quantification RT-PCR of the expression level of IFN-(3 

667 mRNA. The experiment was replicated twice. (G) Vero E6 cells were pretreated with 

668 indicated amount of IFN-J3, infected with wtWIV1, rWIV 1, or rWIV 1 -AX at an MOI 

669 of 0.1 PFU/cell, and posttreated with IFN-(3. Viral replication was analyzed at 24 h 

670 postinfection by plaque assay. The experiment was performed in triplicate and 

671 replicated twice. The differences between selected groups were significant with the p 

672 values less than 0.05: 0.0049 (* bars 4 and 6 in panel A), 0.0008 (** bars 6 and 7 in 

673 panel A), 0.0072 (* bars 4 and 6 in panel B), 0.018 (* bars of rWIVl and rWIVl-AX 

674 in panel F), <0.0001 (*** in panel G). 

675 

676 Fig. 5 Comparison of viral replication efficiency in IFN competent cells between 

677 rWIVl-AX and rWIVl. Calu-3 (A) and HeLa-hACE2 (B) cells were infected with 

678 rWIVl or rWIVl-AX at an MOI of 0.001. Samples were collected at 0, 12,24,36,48, 

679 72, 96, and 120 h postinfection. The viral titers were measured by plaque assay. 

680 
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681 Fig. 6 ORFX protein activates NF-kB. 293T cells were transfected with 100 ng 

682 pNF-KB-Luc, 10 ng pRL-TK, and empty vector, NS1 expressing plasmid, or 

683 escalating (200, 400, 600 ng) of ORFX expressing plasmid. After 24 h, the cells were 

684 treated with TNFa. Dual luciferase activity was determined after 6 h. The results were 

685 expressed as the firefly luciferase activity normalized to that of Renilla luciferase (A). 

686 The relative expression of IL8 mRNA was quantified through qRT-PCR and 

687 normalized to that of GAPDF1 mRNA (B). Differences between selected groups were 

688 significant with p value less than 0.05: < 0.0001 (*** bars 1 and 3 in panel A), 0.0339 

689 (* bars 4 and 7 in panel A), 0.0002 (*** bars 4 and 6 in panel B). n.s., not significant. 

690 The experiments were performed for three times. (C to E) The RNA extracted from 

691 calu-3 cells in Fig. 4 was used for quantification of the expression of IL6 (C), IL8 (D) 

692 and TNFa (E) mRNA. 

693 
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Zeng Page 1 

1 Tables 

2 Table 1 Leader-containing sequences of sgRNAs. 


sgRNA 

ORF(s) 

Leader-containing sequence 3 

CS site 

1 

la/b 

GTAGATCTGTTCTCTAAACGAACTTTAAAATCTGT 

67-72 

2 

S 

GTAGATCTGTTCTCTAAACGAACATGAAATTGTTA 

21,486-21,491 

3 

3aft 

GTAGATCTGTTCTCTAAACGAACTTATGGATTTGT 

25,263-25,268 

4 

E 

GTAGATCTGTTCTCTAAACGAACTTATGTACTCAT 

26,112-26,117 

5 

M 

GTAGATCTGTTCTCTAAACGAACTAACTATTATTA 

26,351-26,356 

6 

6 

GTAGATCTGTTCTCTAAACGAACGCTTTCTTATTA 

26,916-26,921 

7 

X 

GTAGATCTGTTCTCTAAACGAACCACTATGTTACT 

27,272-27,277 

8 

7a/b 

GTAGATCTGTTCTCTAAACGAACATGAAAATTATT 

27,794-27,799 

9 

8 

GTAGATCTGTTCTCTAAACGAACATGAAACTTCTC 

28,300-28,305 

10 

N 

GTAGATCTGTTCTCTAAACGAACAAACTAAAATGT 

28,672-28,677 


3 a Consensus sequence is in bold. 
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